Boundary effects on the scaling of the superfluid density 
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We study numerically the influence of the substrate (boundary conditions) on the finite-size scaling properties 
of the superfluid density p s in superfluid films of thickness H within the XY model employing the Monte 
Carlo method. Our results suggest that the jump p s H/T c at the Kosterlitz-Thouless transition temperature 
T c depends on the boundary conditions. 
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In Ref.§ we used the XY model to investigate 
numerically the finite-size scaling (FSS) properties of 
the superfluid density of superfluid films with respect 
to the film thickness H. The superfluid density p s is 
related to the helicity modulus T of the XY model 
through p s = (m/h) 2 T where m denotes the mass 
of the helium atom[||. We used rather thin films up 
to H = 24 and found that we were able to collapse 
our data for the helicity modulus onto a single curve 
employing the scaling expression 



T(T,H)H eff 
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the effective thickness H, 
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;> where the reduced temperature t 



T/T x - I and 
D (D = 5.8). 

Applying this idea to the experimental data of Rhce 
et al. H we achieved approximate data collapse as 
well (D — 0.145/im). Note that the scaling form 
(Q) reduces to the conventional scaling form in the 
limit H 3> D. In Ref.jjJ we argued that the incre- 
ment D could be understood as an effective scaling 
correction which takes the influence of the boundary 
conditions (BC) into account, i.e. vortex creation on 
a BC dependent length scale I. 

In order to investigate the role of the length scale 
I further we use the XY model on cubic lattices L 3 
with periodic boundary conditions (PBC) along the 
x— and y-directions and Dirichlet boundary condi- 
tions (DBC) (vanishing order parameter) along the 
z-direction and study the FSS behavior of the x- 
component of the helicity modulus with respect to 
L. (In the following we will always refer to the x— 



component of the helicity modulus as the helicity 
modulus.) By studying the FSS behavior of the he- 
licity modulus of the XY model on cubic lattices with 
DBC in the z-direction we preserve the qualitative 
features of films (suppression of the helicity modulus 
due to DBC (cf. also (!]])) but have the advantage of 
using "thicker" systems (up to L = 48). 

The hamiltonian of the XY model on a lattice 
is given by TC — JJ2(ij)^i ' *j where the sum 
runs over nearest neighbors, the pseudospins Si — 
(cos 9 i, sin 9 i) and sit on the lattice sites and J sets 
the energy scale. The helicity modulus T is defined 
as in Refs.Q and the temperature T is measured in 
units of J/ks- DBC are realized on different bound- 
ary inherent length scales as follows. The layers 
z = 1, L are coupled to a spin configuration defined 

by 

s(x,y,z) = (-l)^+Ws(l,l,z), z = 0,L + f, (2) 

with x,y = 1,2, ...,L denoting the integer coordi- 
nates of the lattice sites. The symbol [x/n] means 
the integer part of the number x/n. The integer 
n determines the linear dimension over which the 
spins of the boundary spin configuration are paral- 
lel. Thus, the local magnetization in the boundary 
M = X^ieA i with A an area containing In x In 
boundary spins, vanishes over the length scale 2n, 
i.e. the number n serves as a measure for the length 
scale / over which vortices are created by the bound- 
ary. In our Monte Carlo simulations we used n = 1 
and n — 4 and lattice sizes L = 24, 32, 40, 48. 

In Fig.0 we compare the FSS functions 



T(T, L)L/T = ^{tL 1 ^) {v = 0.6705 @) for PBC, 
n = 1 DBC, and n = 4 DBC along the z-direction, 
respectively. We find ^(sOpBC > ^(^DBC n=l > 
^( X )DBC n=4- This is qualitatively the same be- 
havior as was demonstrated for the film case Q. To 
see the influence of the BC it is instructive to con- 
sider the universal number $(0). We have $(0) = 
0.505(7) (PBC) §, $(0) = 0.189(2) (n=l DBC), 
and $(0) = 0.148(4) (n=4 DBC). Combining these 
results and the results reported in Ref. Q| we can im- 
mediately draw conclusions for the film case. Due to 
the suppression of the hclicity modulus we expect the 
Kosterlitz-Thouless transition temperatures T C (H) 
for films of thickness H in the presence of DBC to 
satisfy the inequality T™ =1 (H) > T c n=4 (ff). Fur- 
thermore the inequality D n —\ < D n —4 should hold 
for the increment D which has to be added to the 
film thickness H to obtain H e ff. Thus, we suppose 
that the larger the length scale I is (or n in our Monte 
Carlo calculations) the larger is the increment D and 
the thicker films have to be used to verify conven- 
tional FSS in experiments. 
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Figure 1: The scaling function T(T,L)L/T for 
PBC, DBC(n=l), and DBC(n=4) along the z- 
direction. 

The dependence of the value $(0) on the BC sug- 
gests that in the case of a film geometry the jump 
T(T C (H),H)H/T C (H) depend on the BC, too. Ta- 
ble | contains the values T(T c (H),H)H e ff/T c (H) 
for films with n = 1 DBC for different film thick- 
nesses and we see that this value seems to saturate 
for H > 16 at about 0.97. Thus, for films with 



H > D we would find T(T C (H), H)H/T C {H) = 0.97 
compared to the value 2/tt expected from renor- 
malization group calculations and found for films 
in the presence of PBC along the z-direction[|| 
and in experiments If we applied n — 4 DBC 
along the z-direction of the film we would expect 
T(T C (H),H)H/T C (H) > 0.97. 
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T(T^(H),H)H eff /T^(H) 


4 


0.565(32) 


8 


0.81(17) 


12 


0.870(65) 


16 


0.974(44) 


20 


0.972(68) 



Table 1: The jump r{T? D (H),H)H eff /T? D (H) 
for different thicknesses H . 

At this point we find it quite interesting to in- 
vestigate the FSS behavior of the superfluid den- 
sity of superfluid films over a wide range of film 
thicknesses and for various substrates so that BC 
are realized which create vortices on different length 
scales. These experiments could directly test our 
ideas of scaling with an effective thickness for 
rather thin films and the dependence of the jump 
T(T C (H),H)H/T C (H) on the BC. It would be espe- 
cially interesting to check these ideas using silicon as 
a substrate and the experimental set-up of the ex- 
periment of Rhee et al.||. 
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